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^\ | In the present paper it is analyzed a very simple example of pseudoextensive system, the tridi- 

mensional system of Linear Coupled Oscillators (OLA Model). The same one constitutes a 
classical tridimensional system of identical interacting particles by means of harmonic oscillators. 
This academic problem possesses a complete analytical solution allowing this way that it can find 
application in modeling some properties of the self-gravitating systems. It is shown that although 
this is a nonextensive system in the usual sense, it can be dealt in the thermodynamic limit with the 
usual Boltzmann-Gibbs' Statistics with an appropriate selection of the representation of the space 
of the integrals of motion. 
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I. INTRODUCTION 



O, 

In our previous works it was established a general methodology to deal with some nonextensive systems. In the 

ref. jjj it was addressed the problem of generalizing the extensive postulates of the traditional Thermodynamics in 

~fs order to extend the application of this theory to the study of some Hamiltonian nonextensive systems. According to 

^- , our proposition, This can be performed taking into account the self-similarity scaling properties of the systems with 

the realization of the thermodynamic limit, the limit of many particles, and analyzing the necessary condition for the 



it was analyzed the most 
ris kind of scaling behavior 



equivalence of the microcanonical ensemble with the generalized canonical one. In the ref. 

familiar self-similarity scaling laws of the systems, the exponential scaling. Systems with tl 

in the thermodynamic limit can be dealt with the usual Boltzmann-Gibbs' Statistics if an appropriate selection of the 

representation of the space of the integrals of movement, 9jv, is taken. That is the reason why we refer those systems 

as pseudoextensive. The extensive systems are just a particular case of the pseudoextensive systems. It is easy to 

show that it is sufficient the presence of an additive kinetic part in the system Hamiltonian for the consideration of 

that system as pseudoextensive. 

The present paper will be devoted to the microcanonical analysis of a very simple nonextensive system that allows 
us a complete analytical study: the tridimensional system of Linear Coupled Oscillators (OLA model). 
This academic model has been used in the modelation of many systems, and although it is very well-known from the 
begining of Mechanics, nowadays it is still considered in the description of some real system, e.L: quantum dots (see for 
example in refs. pH10|). In spite of its simplicity, this model possesses a nonextensive character: it is inhomogeneous, 
the total energy does not scale with the particle number of the system, etc. We will analyze the scaling properties 
of its asymptotic accessible volume in order to precise which is its correspondent asymptotic canonical description in 
the ThL. 



II. THE OLA MODEL 

The Hamiltonian of this system is given by: 
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N N n 

H =E ^p?+ E — (*j - r *) ■ (i) 

3=1 J',fe=l 

This is an example very simple of a self- gravitating small system that allows us a complete analytical study. It will 
develop its microcanonical analysis taking also into account the fallowing integrals of motion: 

N N 

M =y Vj x pj (angular momentum), and Po = /, Pj (lineal momentum). (2) 

3=1 3=1 

In order to eliminate the divergence in the accessible states density of this model due to the unbound movement of 
the system mass center, it is also demand the following additional constrains: 

N N 

P =Y^ Pj = and R =J^ r,- = 0, (3) 

3=1 3=1 

which are consistent with the energy and angular momentum conservation. It is convenient to work with dimensionless 
variables, so that, it is introduced the following units: 

[E]=Huj, [M] = ?i, [R]=( V and [P] = (mfiw)* • ( 4 ) 

\rauj J 

The calculations will be facilitated using the vectorial convention for the R 3Ar vectorial space, which appears in the 
appendices. It is easy to see that the equations Eq.([j]), Eq.(g) and Eq.(||) are rewritten as: 

H =\p 2 + ^(b 2 R 2 -(B-R) 2 Y M = R®P, (5) 

Ro = B ■ R = and P = B ■ P = , (6) 

where R,Pe R 3N are the extended system coordinates and linear momentum that together conform the TV-body 
phase space of the system, and B E R N with B = (1, 1, . . . 1). The solution of this problem can be easily found 
through of the generating functional of the distribution: 

1 f d 3N Rd 3N P 
K( X ,v;I,N) = -t- e X p[i( X -R + vP)]6[I-I N (R,P)], (7) 

h LU J (27r) 

where / = (E, M, R , P ), \ and v are 3-/V-dimensional vectors belonging to R 3Ar . It is convenient to use the Fourier's 
representation of the delta function: 

f + °° dk 
d [x - x) = / — exp [z(x - x)\ , (8) 

J-oo 27T 

where z = e + ik with e > 0, to rewrite the Eq.(^) as follow: 

1 f <i 10 k 
K( X ,v;I,N) = - r - / — - w exp (z • /) H (z; x , v) , (9) 

a uo J (2-k) 

where k = {k\,\z-2,ip,iri), and z = (zi, Z2,ip,ir)), with z\ = /3 + iki, z 2 = 7 + ik2. The function H (z; Xi v ) is defined by: 

f d 3N Rd 3N P 
K(z; X ,v)= / exp{-[z- I N (R,P) - i X - R- iv ■ P}} . (10) 

J (27r) 

The above integral is very easy to calculate because it can be reduced to gaussian integrals. The calculations yield: 

+00+00+00+00 
K(x,v,E,M,N)=^J_ 7 f f f f * A -p(^)exp(z 2 .M) 



(2TryB^h'u J J J J (2^) 4 (Bz 1 r- 1 (B*z*-3&y 



— co — 00 — CO — 00 



/ B*v 2 -{B-v) z \ / B 2 z 2 A 2 - (Bz 2 ■ A) 2 \ / B 2 (z 2 x A) 2 - (£z 2 ® A) 2 \ 
X 6XP ( 2^ J CXP ( 23^ J ° XP ( 2£ 2 z 2 {WzJ - z 2 ) Zl ) (11) 

where: 

A = X +^ (12) 

The above expression must be used in the obtaining of the physical observables. The accessible states density is 
expressed as: 

+OO + OC + OC + OC 

q(e,m,n)= 3 1 7 I I I I dkld y ^f)^ 2 '^) (13) 

{2nfB^h 7 oj J J J J (2tt) 4 (Bzxf-UBHl-zlf- 1 



— oo — oo — oo — oo 



whose integration yields: 



The function Hn (x) is given by: 

-Hjv (a;) = (1 — %) Pn-2 (1 — x) a {x — 1) , where cr (x) = < n . f ^ n is the Heaviside function, (15) 

and Pm_2(z) =\^ a n z n , with a„ = A— — --r- rr- t- (16) 

w ^ (iV-2-n)!(2JV-4 + n)!(7V-l)! 2™n! v ' 

The accessible volume is obtained multiplying the state density by an appropriate element volume constant SI. A 
reasonable choice for the tridimensional OLA model is: 

SI = h 7 uj, (17) 

and therefore: 



W(E,M,N) = — * M3JV-7# (JL\ 
{ ' (2nB) 7 \BM) 



(2ttB)' 
The integrals of motion space, ^sn, in the H-ie,'m,n) representation is the R. cone for a given N: 

E 2 - AM 2 > 0. (18) 

In order to access to the scaling laws, it must be obtained the asymptotic dependency of the polynomials coefficients 
in the Eq.(|16[). It is easy to show the following asymptotic behavior: 

1 (-l) n A" 

hm a„ ~ —^T7 — ^tt : — ■ (19) 

jv-oo 2 2N N 3N 4 n n! 

Therefore, the asymptotic accessible volume is given by the expression: 

Wasym(E,M,N)<X \— ^-g- ) Qn(^^), (20) 

V.2(JV)V v E J 

where the function Qn (x) is expressed as: 

Q N {x) = (l-x) 2N exp(N^-— -J , with x € [0, 1] . (21) 

It is very easy to see that the tridimensional OLA model possesses an exponential self- similarity scaling laws [0,0: 



N -> N {a) = aN 

E -> E (a) = al E } => W asym (a) = T [W a , a] , (22) 

M -> M (a) = aM 

where the functional T \W Q , a] is given by: 

T [W , a] = exp [a In W ] . (23) 

where a is the scaling parameter and W is given by: 

W = W asym {E 0l M ,N ). (24) 

Thus, the tridimensional OLA model is not extensive, but pseudoextensive [|2|. In the present case, the generalized 
Boltzmann's Principle W adopts its ordinary form: 



S B =lnW, (25) 



his celebrated gravestone epitaph in Vienna. To this entropy form corresponds the usual Shannon-Boltzmann-Gibbs' 
extensive entropy: 

Ssbg = - ^2 Pk lnpfc ' ( 26 ) 

k 

and therefore, in the thermodynamic limit, the microcanonical description could be substituted equivalently by the 
usual Boltzmann-Gibbs' Distribution with a appropriate selection of the representation of the integrals of motion 
space. This can be performed changing the representation from TZ(e.m.,n) to 1Z(s,m.,n): where: 

£ = -=E. (27) 



In the thermodynamic limit the microcanonical description of this system should be substituted equivalently by 
the canonical one. In this case the Boltzmann-Gibbs' Distribution is: 

uj c (X; [3, 7 , N) = 7(R l exp [-[3£ N (X) - 7 • M N (X)} . (28) 

Through the Laplace's Transformation: 

Z ((3, 7 , N) = J exp {-[3E - 7 • M) W (£ , M, N) ^M (29) 

it is easy to show that in the thermodynamic limit both descriptions are equivalent, that is, the Legendre's Transfor- 
mation between the thermodynamic potentials of the ensembles is valid: 

P(/?,7,A0~/?£+7-M-S s (£,M,A0, (30) 

where P (f3, 7, N) is the Planck's Potential: 

P(/?,7,A0 = -lnZ(/3, 7 ,A0, (31) 

when exists a unique sharp maximum in the integral argument of the Eq.(E9). In this maximum, (£ m ,M m ), it will 
be satisfied the following conditions: 

13 = -^Sb (£ m ,M m ,N) , 7 =g^S B (£ m ,M m ,N) , (32) 

and all the eingenvalues of the curvature tensor Mm : 

K^ = d^d v S B (£ m ,M m ,N), (33) 

are negatives [in the Eq.(p3) d^ = gj^-, where 7 M = (f ,M)]. Similarly, the generalized Duhem-Gibbs' relation is also 
valid in the thermodynamic limit: 



S B {£, M,N)={3£ + ~/- M+fiN, (34) 

where fi is the chemical potential: 

n(0,"/) = -^S B (£ m ,M m ,N). (35) 

From the Duhem-Gibbs' relation is easily deduced the following relationship: 

P(f3, 1 ,N) = -n(/3, 1 )N. (36) 

The Boltzmann's entropy of the tridimensional OLA model in the asymptotic region is given by: 

and therefore, the canonical parameters (3 and 7, and the chemical potential [i, are given by: 

l Ar 12£ 2 -3£Af-M 2 1 9£-M M S B (£,M,N) . .... 

= 4^ g 2 (£ - M) ' 7 = ~4 N £{£—m M' M = iV 3 ' (38) 

To validate the previous results it must be demanded the concavity of the Boltzmann's entropy, that is to say, the 
non-negativity of the 4x4 curvature tensor: 

(1 at 6£ 3 -3£ 2 M+M 3 1 Ar 9£ 2 -2£M+M 2 M 3x i \ 

2 JV £ 3 (£-M) 2 i £ 2 (£-M)* M \ ,oq\ 

1 at 9£ 2 -2£M+M 2 M 3x i » f 2 M 3x i Mi x3 1 at 9£-M /t \ I V / 

4 £ 2 (£-M) 2 M (£-M) 2 M M 4 iv £(£-M)M l*-*»V3x3 ) 

where Mi X 3 = (M3 X i) = (M x , M y , M z ), and Ij r is the 3x3 transverse unitary matrix: 

The determinant of the curvature tensor is given by: 

detX = -L^ 15g2 + 18g 4 M ' M2 0* - M) 2 > 0. (41) 

256 (£-Mf£H4 2 V ' y ' 

From the above result is derived that all the accessible space is appropriately described by the canonical ensemble 
in the 7£(£,m,jv) representation: in the Laplace's Transformation, Eq.(E3), there is only one sharp peak. Thus, in the 
thermodynamic limit both descriptions are identical: all the system accessible space can be deal with the canonical 
description. This result facilitates so much the analysis of this system. In the canonical description, the Planck's 
potential is obtained from the Eq.([ll|) replacing Bz\ — ► (3 and Z2 — > 7: 

P(/3, 7 ,iV) = ^ln[/3(/3 2 - 7 2 )]. (42) 

The generating functional in the canonical description: 

r d 3N Od 3N P 
G c ( X , v- /3, 7, N) = / * N lo c (Q, P- [3, 7, N) exp (i X ■ Q + iv ■ P) , (43) 

J (27T) 



is also obtained from the Eq.(ll) 



/ B^ - (B ■ vf \ I B^\^-{B 1 -\) I \ I J3 2 ( 7 xA) i -(i37^A) i ^ 
G C (X,V,^,N) = exp ^ — J exp ^ ^-^ j exp ^ ^fW^f) * j " 

(44) 

It is very interesting to known the particle distribution in the space. This observable is easily obtained from 
the generating functional as follow: setting zero all the components of the vectors v and leaving only one vectorial 
component in the vector Xi X = (li 0; ■■■■7 0); an d performing the following integration: 



+oo-|-oo+oo 

dqid 2 q t . ( B 2 - 1 „\ / B 2 - 1 



/' ( '■ I = V / / / :r g^T ex P i-inqi) exp (-ir t • q t ) exp \-—^Zq*j exp ^-^-^— ^y/3q 2 ) , 1 45 ) 



— oo — oo — oo 



where qi and q* are the longitudinal and the transverse components of the vector q with respect to the vector 7. The 
integration yields: 

W / r 2 2 \ 

'M-vpssT* (-£-£)■ (46) 

showing that the particle distribution possesses a gaussian shape where the parameters 07 and at are given by: 

B 2 — 1 B 2 — 1 

Let the parameter s be introduced as follow: 

parametrizing a family of ellipsoids of revolution which characterizing the surfaces of equal density: 

N f s 2 \ 

P0) = ^^=exp -— , (49) 

where a is defined by the relation: 



o=\Joio 2 - (50) 

The eccentricity of ellipsoids family is constant and given by: 



*=f-i-h-E- (51) 

As it can be seen, the rotation deforms the spherical shape of the distribution. The characteristic size of the 
distribution is the given by: 

lc = V^- (52) 

From the Eq.(f47|) it is deduced that the characteristic size decreases with the increasing of the particle number N: 

l c ~l/<fN. (53) 

This kind of scaling law is not typical in the extensive systems, which are usually scaled proportional to the systems 
size. In the FIG.l. it is shown different profiles for the particle distributions. 

III. CONCLUSIONS 

In the present paper we applied the results of our previous works on study of a very simple model of self-gravitating 
small system, the tridimensional OLA model. The microcanonical analysis showed that the OLA model is a pseudoex- 
tensive system pj, since it possesses a exponential self-similarity scaling laws in the thermodynamic limit. Systems 
with this kind of scaling laws could be dealt with the usual Boltzmann-Gibbs' statistics if an appropriate represen- 
tation of the system integrals of motion space is selected. Although the model presented before is very simple, the 
analysis shown that it is necessary the consideration of self-similarity scaling postulates for the analysis of this system, 
since the same is not extensive. 



IV. APPENDIX: THE VECTORIAL CONVENTION. 

Let R 3N be the 3./V-dimensional vectorial space, which is the external product of the spaces R 3 , R N . The bases of 
R 3JV are represented by means of the product of the bases of the respective spaces, R 3 and R^: 

E M = t.~E k . (54) 

It will be only considered those transformations preserving this property, that is, those transformations which acts 
in the spaces R 3 and R N separately. The vectors of this space can be represented as: 

N 3 N 3 

v =J2Y1 x k, s ~?s~$k =Y. ^ k ^ k =Y1 ^»~?» where ^ eR3 ;^ eR " • (55) 

fe=lS=l fe=l 8=1 

From here and after, it will be adopted the tensorial summation convention. Using the above representation, it can 
be defined the following external operations: 

Let be the vectors: o> G R 3 ;^G R^V^i,^,^ € R 3JV - 

Partial and Total Inner Operations: 

R 3 x R 3N -» R N : -ct ■ V = (1? • lf k ) ~^ k , (56) 

R*xR 3 ^R 3 :l.y=(l. V s ) "e%, (57) 

R 3 x K 3N - R N : V ■ W = (vt ■ vf k ) - (ft ■ W s ) . (58) 

Partial Vectorial Product: 

R 3 x R 3N -* K 3N :lt xt = (ltxlf k )^ k . (59) 

Vectorial-Scalar Product: 
R 3JV $ R 3N -+ R 3 : T? O M? = (^ x uf fc ) . (60) 

From the previous definitions is derived the following identity: 

-ct-\?®W =(lt xV^j-f . (61) 
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